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1. INTRODUCTION 
Importance 

The United States possesses son« of the most productive agricultural 
lands In the world. Indeed, the millions of hectares of land In agricul- 
tural production In the U.S* are an Important natural resource. Mainten- 
ance of the quality and quantity of agricultural land, as a ^'enewable 
resource, Is Important If we are to maintain or improve our future balance 
of payments picture. In 1980, the United States exported over $40 billion 
worth of agricultural products. Since overall U.S. exports In 1980 totaled 
just over $200 billion, the portion which agriculture contributes Is sig- 
nificant, It Is even more significant when one realizes that agriculture 
consistently produces trade surpluses. In 1900 we exported over $21 
billion more agfl cultural pr-oducts than we imported. Non -agricultural 
export categories showed a trade deficit of $44.8 billion. So with our 
agricultural surplus our total foreign trade deficit for 198(> narrowed 
to $22.6 billion. 

United States agriculture today is beset by a growing number of 
problems ranging from expansion of urban areas to dwindling water supplies. 
In 1972, the U.S, Departnent of Agriculture reported that during the 1960 ‘s 
some 295,000 hectares of land were urbanized each year; recreational land 
uses increased annually by some 410,000 hectares while transportation re- 
lated land uses expanded at an annual rate of 55,000 hectares (USDA, 1972). 
Recent figures indicate this trend has continued and has Intensified, 

In addition to the problems of land conversion, U.S. agriculture Is 
faced with water problems as well. Conservation of water resources Is not 
just a U.S. problem; It Is a truly world-wide problem. Competition for 
water is increasing everywhere between agricultural, urban, Industrial, and 
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t^creational users. In many areas of the globe demands for water are ap- 
proaching or, in some cases, exceeding supplies. The major user of water 
in the world is irrigated agriculture. 

In the United States, agricultural irrigation accounts for over 80 
percent of the water consumed (GAO report RED-76-116, 1976). Ninety 
percent of irrigated agriculture in the United States is in 17 western 
states where over 20 million hectares (50 million acres) are irrigated 
(Bureau of Reclamation, 1976). From declining groundwater tables in high 
plains aquifers to the need for new transportation, delivery, and drainage 
systems in Sun Belt states, water availability is an increasing problem. 
Efficient use of available water resources is becoming increasingly im- 
portant. This is particularly true in southwestern Sun Belt states where 
commercial irrigation agriculture is practiced on a large scale. In the 
future, as population in these states continues to rise .and the costs of 
new irrigation projects also rise while construction of approved projects 
are slowed by the need to satisfy environmental legislation, competition 
for available water supplies will become pven more intense. 

Because of the increasingly competitive demands placed on available 
water resources, irrigation agriculture must be made more efficient. 
Irrigation agriculture must move toward more optimum use of its share of 
our increasingly short water resource. The need for work in this area 
is even more important if we believe the 1976 Government Accounting Office 
(GAO) report which stated that "over half of the water delivered to farms 
is wasted through overwatering which can limit crop production, increase 
farming costs, and contribute to water pollution" (U.S. Gov't. 1976). 

This same GAO study cited a number of instances in which additional lands 
could have been put into production if irrigation water had not been used 
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excessively on lands already under production. It is our belief based 
on the research presented herein, that remote sensing can make a contri- 
bution toward more efficient use of available water resources and for 
services provided In monitoring water allocation. 

Background 

The U.S. Bureau of Reclamation (USBR) Is responsible for planning, 
constructing, operating, and maintaining facilities for storing and 
delivering water supplies for Irrigation in the western U.S., and thus 
has an Interest In promoting efficient irrigation management. In 1976, 

USBR delivered over 26 million acre-feet of water to 3.6 million hectares 
(20 percent of the irrigated land in the western U.S.). 

In 1970, USBR began developing the Irrigation Management Services 
(IMS) to assist Irrigators toward mom efficient use of water through 
Irrigation scheduling. IMS provides two levels of scheduling assistance 
to the farmer - the “Irrigation Guide", and the “Field Irrigation Schedule". 
The Irrigation Guide gives irrigation Intervals, based on dally estimated 
evapctransplratlon (ET) rates and average water-holding capacity of major 
soils, for principal crops in an area. The Field Irrigation Schedule 
provides the farmer with up-to-date soil moisture status for each of his 
fields In the program. The schedule provides recommended Irrigation dates 
and amounts of water to be applied. One field technician is currently 
required for every 800 to 1600 hectares served by this program. 

In 1977, 63,516 hectares (156,884 acres) were served at the field 
level and approximately 35,000 hectares (87,000 acres) served at the guide 
level. When combined with ancillary services provided by IMS, a total 
of 142,000 hectares (350,000 acres) were Involved in Irrigation scheduling. 
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The field-level service Is based on a procedure that uses dally 
estlniates of ETi soil moisture, effective precipitation, and applied 
Irrigation water to compute dally soil-water depletion (Jensen, 1968). 

By projecting the ET rate Into the future, the date at which soil -water 
depletion reaches Its maximum allowable value, Is determined. 

Actual dally ET for a given crop Is estimated from potential ET 
(calculated using meteorological data and a modification of the Penman 
(1963) combination equation) and a crop coefficient Kc, defined as 

Kc - ET/ETp 

where Et Is actual evapotranspi ration and ETp Is potential ET, The crop 
coefficient Is Influenced by percent ground cover, leaf area Index (LAI), 
growth stage, and soil moisture (Ritchie and Burnett, 1971 j Kanemasu and 
Powers, 1974), all of which have potential for measurement using remote 
senslng- 

Crop ooefficlent curves used by IMS are developed from experimental 
data and represent ET of crops In a full and healthy population. Field 
evaluations of crop growth are used to estimate crop coefficients from 
experimental curves such as those shown In Fig. 1. ET rates estimated 
using experimental curves require an adjustment because crops growing 
under experimental conditions generally are not representative of crops 
growing under normal Irrigated conditions. Field measurements of soil 
moisture are currently used to adjust Kc for the effects of soil moisture 
(Buchheim and Ploss, 1977), At present, only tedious ground sampling 
procedures are available to acquire these necessary inputs. 

Spectral Information acquired on a repetitive basis can potentially 
provide Inputs into the IMS program. Visible and reflective Infrared 
wavelengths can be used to evaluate vegetation parameters (LAI, percent 
/jover, etc,) which affect Kc (Kanemasu et al . , 1974; Kanemasu et al.. 
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1977; Wiegand et al., 1977; Tucker, 1979), Thermal Infrared and micron- 
wave wavelengths can be used to evaluate near surface soil moisture as 
surniiarized by the NASA/USDA Soil Moisture Workshop (Heilman et al., 1978). 
All of the remote sensing approaches have certain rudimentary details 
available in the literature, but for an operational system for all crop 
growth/soil moisture conditions, a detailed evaluation must be conducted. 

The Remote Sensing Section of USSR has a computer-interactive image 
analysis system that can be integrated inii* an operational system of 
irrigation scheduling. In addition, IMS has an already established in- 
formation dissemination network to rapidly transfer information to the in- 
dividual farmer. Both these aspects are critical to the actual implemen- 
tation and use of any developed procedures. 

A grant was established 1 April 1980 with the long term objective 
of testing the utility of thermal and reflective data for predicting crop 
coefficients and soil moisture as inputs into the U.S. Bureau of Recla- 
mation irrigation scheduling program, The grant proposal was in response 
to RTOP 677-'^2-23 entitled Irrigation Scheduling. 

A ground study utilizing hand-held radiometers was to be conducted 
at the Navajo Indian Irrigation Project near Farmington, New Mexico in 
conjunction with ongoing USBR measurements including lysimetric measure- 
ments of evapotranspi ration. The Remote Sensing Institute (RSI) was to 
collect visible, near IR and thermal IR measurements to be analyzed and 
evaluated in terms of the ground measurements which included percent crop 
canopy cover obtained by RSI. Results were to be used to recommend future 
action regarding the use of satellite data in irrigation management. 






2, RESEARCH PROGRAM 
Inth)duct1on 


RSI conducted a field measurement program during the summer of 1980. 
Measurements were made with two visible/near IR hand-held radiometers 
and a hand-held thermal radlrmster, Soil moisture and lysimetric measure- 
ments were collected by USBR personnel and forwarded to RSI along with 
meteorological data. 

Field Measurements 

The study was conducted on the Navajo Indian Irrigation Project near 
Farmington, New Mexico, In conjunction with a Bureau of Reclamation water 
use study. The test site consisted of four alfalfa plots with lysimeters 
that were irrigated with a line souna sprinkler. The line source linearly 
decreases the amount of water applied as the distance from the source 
Increases, With the configuration of the system, plot 1 was over-irrigated, 
plot 2 received optimum applications, and plots 3 and 4 were under-irrigated. 

Radiances from each plot were measured with Mark II three-band and 
EXOTECH four-band radiometers in the m1d-morn1ng (1000-1100 MDT) and early 
afternoon (1300-1400 MDT). Radiances were converted to reflectances using 
radiance measurements of reference panels coated with 3-M v^hlte Nextel 
paint and BaSo^. Panel radiances were measured before and after each plot 
measurement. Radiometric temperatures of the plots were measured with a 
PRT-5 within the same time frame. Radiances of bare soil were also measured. 

Surface soil water contents (0^4 cm) were determined gravl metrically 
on samples collected at the same time as the spectral measurements. Per* 
cent canopy cover was determined using 35mm color slides (photographed 
from a vertical position) projected on a random dot grid. All measure- 
ments began after the second cutting of alfalfa in eey'ly July 1980 and 
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continued to maximum canopy cover. 

USBR personnel coll ectet profile soil moisture (neutron attenuation 
method) and lysimetric measurements of actual evapotranspi ration. The 
collected meteorolooical data included maximum and minimum air and dew 
point temperatures, solar radiation, wind speed, and pan evaporation. 

USBR calculated potential ET using the Penman equation which was cali« 
brated for the climatic conditions at Farmington. 

Preliminary Results 

wtm 

The scheduling procedure used by the Irrigation Management Service 
(IMS) of USBR requires estimates of the crop coefficient which is the ratio 
of actual to potential ET. We evaluated the relationship between Mark II 
spectral measurements and the crop coefficient. 

Best relationships were obtained with the perpendicular vegetation 
index (PVI) of Richardson and Wiegand (1977) using channels 1 (0.63-0.69|jm) 
and 2 (0.76-0, 90pm). The soil background line was 
channel 1 * 1,03 channel 2 - 0,05 
with an r^ of u.89. 

Figures 2, 3 and 4 show the relationships between PVI, the percent 
cover, and the crop coefficient. A PVI based on channel 2 and 3 (1 .55-1 ,75pm) 
did not improve the results. 

These results demonstrate the potential for using spectral measure- 
ments to estimate the crop coefficient. Inclusion of a liquid water 
absorption channel does not appear to offer a significant benefit. 

The Appendix is a paper submitted for publication which resulted from 
this project. The paper, "Evaluating the Crop Coefficient Using Spectral 
Reflectance," submitted to the Agronomy Journal , is the culmination of the 


Mark II Radiometer data analysis activities. 
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Figure 2. Relationship of the crop coefficient K (ET/ETp) to percent cover for 
1980 alfalfa plots. 
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Figure 4. The relationship of PVI and K for the alfalfa plots. 
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3. FOLLOW-UP TO THE RESEARCH PROGRAM 
The Navajo Indian Ir^rlgation Project 

The Navajo Indian Prx)ject near Farmington, New Mexico, which encom- 
passes the study area, Is not presently an operational program— It is 
a consumptive water use study. The long term plan at the location is 
for 105,000 acres to eventually be under irrigation management. Two 
aspects being considered in this large scale Irrigation management program 
are: 

1) water diversion requirements within the main water 
conveyance system, and 

2) on-farm scheduling for water. 

Although on-farm scheduling is being discussed, only water diversion re- 
quirements are presently being addressed. There is one turn-out from 
the main system per farm. The water diversion program now assumes that 
after water is delivered to the turn-out in response to a water order from 
the farmer to the irrigation district that those farm fields were irrigated. 
It further assumes that another official water order will come a certain 
number of days later. There is no more detailed interaction between 
the district and the farmer than the water order. 

On-farm scheduling would require a much more extensive interaction. 

The water district, through computerized modeling of on-farm water usage 
could more accurately determine its scheduling, but would require more 
information about the individual farms and their usage. Computerized 
modeling would require many inputs— crop type and acreage, soil moisture 
and others, Present discussions with the Bureau of Indian Affairs and 
the Navajo Indian Irrigation Project personnel about on-farm scheduling 
do not include modeling and concern using neutron probes to obtain soil 
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moisture data. 

The possibilities for interfacing remotely sensed crop coefficient 
data into an on-farm scheduling service were discussed with Jerry Buchheim 
and Bob Hanson, Bureau of Reclamation, Denvei". The results of the research 
effort were presented, illustrating the potential for estimating crop 
coefficient from Landsat-type spectral data. The Bureau's information 
needs toward development of on-farm scheduling include crop inventory 
techniques and consumptive water use monitoring through the growing season. 
Proposed application projects utilizing remote sensing must be submitted 
through the cooperating irrigation district and require a monetary con- 
tribution by the district in order to obtain Bureau approval. The idea 
of a demonstration project was appealing and it was agreed that future 
cooperation toward such a project was desirable. 

Truck-mounted Radiometer System 

The remotely sensed data utilized in the crop coefficient estimation 
portion of the grant effort were acquired with hand-held radiometers. 

RSI was given the opportunity to upgrade considerably it's capabilities 
for multispectral measurements on the ground and has done so with a little 
help from this grant although it was not one of the stated objectives 
of the project. RSI was designated by NASA as one of a few institutions 
which best could take advantage of the new LARS/Barnes modular multi - 
band radiometer in on-going NASA-sponsored projects. 

RSI purchased a 3/4 ton pickup. LARS/Purdue manufactured a boom 
and mounting platform and a separate reference panel platform for the 
RSI truck under one of their NASA contracts through the Johnson Space 
Center. RSI utilized some travel funds from this NASA Grant to transport 
the truck to and from Purdue for the installation of the boom and platforms. 
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The LARS/Barnes Radiometer has not been received at RSI yet, although 
the other elements of the data system have. Therefore, no actual measure- 
ments are available with which to illustrate this new capability. The 
system will be utilized in South Dakota during the 1982 growing season 
as the primary data acquisition system on one or more experiments con- 
cerning remote sensing of cultivated crop parameters, 

4. A LOOK TO the FUTURE 

An outgrowth of the review of source material and analyses of the 
data acquired is the realization of the potential for developing a remote 
sensing assisted Irrigation Scheduling Information System (ISIS). Such 
a system would contain data gathered from a variety of sources depicting 
information for regions of differing sizes all registered to a common 
cartographic base with a uniform scale (See figure 5). Different data 
planes could be input to the system at regular or irregular intervals. 

ISIS is conceived of as a flexible system Which assumes a high degree of 
user interaction and indeed cannot adequately function without user 
participation. That is, any remote sensing assisted irrigation scheduling 
information system can only work if individual irrigators (farmers) fully 
cooperate. Individual farmers therefore must pain an understanding of 
and an appreciation for the fact that the data they supply must be accurate 
and tijTiely in order to insure that the information produced by the system 
is useful to them in terms of optimum scheduling. ISIS should not 
necessarily be considered an isolated system tut could, and probably 
should, be considered more as an integral part of a larger farm management 
information system such as the Agnet system developed at the University 
of Nebraska. 

Paralleling the Agnet systems concept it is envisioned that ISIS 
be comprised of a central computing facility with terminal links to 






^ J ta t *U 


iW 

< »5 i® • 

lill: 




tU 

> 1 * 
ri 
1 4 

%<4 »! 

35:33 


3 U) 

■e 

m 

iH 


'H1“:^S *(S.5 
"rl 

<^<^1 -iNlS 


?/2ii 

u.' (1 
^ 2 


nmmM 


QlRIGlNM. PASU w 
OF POOR QUAUW 


9 

u 

'll 
^ 9 3 5 
v3ip rS 

5 S 5 1 1 

S“e !|-|4 

Q * ig S 5 ^ 

;j!S».a If 


a. 0 ?• « 

i.M4 

fiJgS^^bSoS^c f 3 

Is^-P 

S4- 

”4||lsl|l|||| 

SuSSp 5<03««““'^ 

i 




farm agents' offices and eventually to individual cooperating Irrigated 
farmsteads, Table 1 shows the Information vhich would be held in each 
separate location and that which would be shared by both computers, 
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TABLE 1 


INFORmTION RESIDENT ON ISIS COMPUTERS 


Farmer Computer Shared information 


Central Computer 


Field Boundaries 
Soils 

Farm Botiidaries 
Location of Tile 
Lines 


Crop Type 
Planting Oates 
Initial Soil moisture} or 
Date of Pre-Irrigation 
Irrigation Canals 
Gates 

Drainage Ditches 
Irrigation District Boundaries 
Other Political Administrative 
Boundaries 
Stage of Growth 
Date of Last Irrigation 
AmouR''; of Last Irrigation 


Leaf Area Index 
% Vegetation Cover 
Meteorological Conditions 
Water Budget Calculation 
Crop Condition Calculation 


By distributing Information storage In this tiwnner and implementing coded 
access procedures! security of Individual farmer Information could be main 
talned. Such a network also means that the farm computers could be small 
and the access time and costs mlnimllted by having the Image processing 
and modelling routines which require more computing power and frequent 
updates residing on a larger central computer (possibly in a regional 
Bureau of Reclamation office). 

The ISIS concept, as presented herein, is not Intended to replace 
the current operational Bureau of Reclamation irrigation scheduling pro- 
cedures described 1n the background section of this report. We do, how- 
ever believe that implen«ntat1on of a cooperative research program with 
the Bureau to prove the effectiveness of an ISIS type concept could pro- 
vide significant benefits. We believe that an Irrigation scheduling In- 
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formation system could facilitate the work of Bureau field-level personnel, 
But the real beneficiaries of such a system would be the farmers (irri- 
gators), By tapping into the geobased information potentially available 
from an ISIS type system, the individual farmer could fine tune (in 
essence, customize) his irrigation procedures for each of his fields, 

Such a system can improve irrigation efficiency and potentially save 
water while maintaining and hopefully increasing yield through improved 
management and/or bringing new lands into production. 
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"Evaluating the Crop Coefficient Using Spectral Reflectance" 
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ABSTRACT 


A field study was conducted In four differentially Irrigated 
plots of alfalfa ( Medicago sativa 1.,) planted In Shiprock sandy loam 
to assess spectral reflectance for estimating the crop coefficient 
(Kc)» the ratio of actual to potential evapotranspi ration (ET). A 
bidirectional reflectance factor was measured using a three-channel 
(0.63 to 0.69 ym, 0.76 to 0.90 ym» and 1.55 to 1.75 hand-held 
radiometer, and was used to calculate a perpendicular vegetation Index 
(PVI). Actual ET was determined by the water balance method In non- 
weighing lysimeters, and potential ET was calculated using Penman’s 

t 

equation. 

Significant linear relationships were found between PVI, and per- 
cent cover (r^ * 0.911), and between Kc and percent cover (r^ * 
0.815). In addition, the position of the PVI Intersection on the soil 
background line changed as a result of soil moisture increases follow- 
ing Irrigation, even at high percent cover. Thus, once experimental 
relationships between Kc and crop growth are established, a mean 
Kc can be determined from spectral estimates of stage of development 
and the soil background component of PVI can be usi^d to adjust the 
mean Kc for Increased evaporation following Irrigation because the 
ratio of actual to potential evapotranspi ration will approach 1 when 
the soil surface Is wet. 
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Most Irrigation scheduling models have an evapotransplratlon (ET) 
component to calculate soil water depletion# A common method of esti- 
mating ET Is to relate it to potential evapotransplratlon (ETp) us- 
ing a dimensionless crop coefficient Ko defined by the equation 

Kc - ET/ETp. Cl] 

The scheduling models may differ In the method used to determine Kq 
olid ETp • 

The crop coefficient Increases from a low value at emergence to a 
maximum s’* full cover. However, significant fluctuations In the ET/ 
ETp ratios can occur due to Increased evaporation from the soil sur- 
face following rainfall or Irrigation. Kq will approach 1 when the 
soil surface Is wet, regardless of crop stage of development. In 
addition, Kq can be Increased by advection, or decreased when soil 
moisture becomes limiting. 

Experimental crop coefficients have been determined for a number 
of different crops, and are usually e>. "essed as a function of growth 
parameters such as percent cover, leaf area index, or growth stage 
(Denmead and Shaw, 1959; Fritschen and Shaw, 1961; Grimes et al., 
1969; Jensen, 1968; Jensen et al., 1970; Kanemasu and Powers, 1974). 
Use of the experimental data for scheduling Irrigation on a f1eld-by- 
fleld basis Is difficult since fields of the same crop species are 
seldom at the same stage of development. In addition, the experimen- 
tal relationships do not account for Increased evaporation following 
irrigation. 
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Because Kc Is closely related to crop growth, It should be 
possible to use spectral reflectance to derive a crop growth parameter 
for defining a crop coefficient curve. Development of remote sensing 
techniques to estimate K<; would be particularly beneficial for large 
Irrigation projects where manpower limitations preclude frequent field 
observations necessary to accurately determine Kc and ET. 

We conducted a field Investigation in differentially Irrigated 
plots of alfalfa ( Medi cago s^:iva L.) to assess the use of spectral 
reflectance for estimating the crop coefficient, and for adjusting the 
crop coefficient to account for Increased soil evaporation following 
Irrigation or rainfall. 
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Experiments were conducted on the Navajo Indian Irrigation Pro- j 
ject (NIIP) near Farmington, New Mexico. Average annual rainfall of 
the area is 176.3 mm. 

The research site consisted of four 1.8 x 1.8 x 1.8 m lysimeters 
in a field of Shiprock sandy loam (ccarse-loamy , mixed mesic Typic 
Haplargid). Lysimeters were installed so that tops of the walls 
were 15 cm below the soil surface. Each lysimeter contained two neu- 
tron probe access tubes. Lysimeters were drained with perforated pipe 
and ceramic cup suction systems. Alfalfa (Medi cago sativa L.) was 
planted in the 32-ha field in late July 1979. 

The field was irrigated with a line source sprinkler system 
(Hanks et al., 1976) that produced a continuously decreasing irriga- 
tion pattern at right angles outward from the line source. The 
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lysimeters (designated plots 1 through 4) were U5, 4.5, 7.5, and 10.5 
m, respectively, from the line source. The plots were irrigated twice 
during the study (Table 1). 

Spectral radiances were measured with a hand-held 15 degree 
field-of-view radiometer (Tucker et al., 1981) in three wavelength 
bands (0.63 to 0.69 ym, 0.76 to 0.90 ym, and 1,55 to 1.75 >jm) at a 
height of 1.5 m above the canopy. Data were acquired between 1000 and 
1200 MOT following the first alfalfa cutting (7 dulyi 1980). On each 
date, six measuremer-s were made per plot and averaged for analyses. 
Bare soil radiances were also measured in a similar manner on a plot 
adjacent to the alfalfa. 

A panel coated with highly reflecting barium sulfate was used for 
determining the bidirectional reflectance factor (Silva, 1978), 
defined as 

R.(X) 

BRF(x) " ITTxy * Pr^^^ 

where BRF(x) is the bidirectional reflectance factor at a specific 
wavelength interval x, Rs(x) is the radiance of the scene, Rp(x) 
is the radiance of the reference panel, and pp(x) is the reflectance 
of the reference panel at x* Radiance of the reference panel was mea- 
sured immediately before radiance of each plot was measured. 

The BRF*s were used to calculate a perpendicular vegetation index 
(PVI) (Richardson and Wiegand, 1977), defined as 
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PVI - [(Rgl - RPi)^t (Rg2 - Rp2)^]^/^ [4] 

Where Rpl and Rp2 are composite BRF's (vegetation plus soil back- 
ground) measured in radiometer channels I (0.63 to 0.69 pm) and 2 
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(0.76 to 0.90 Jin), respectively, and Rgl and Rg2 are soil background 

\ 

BRF's corresponding to the vegetation points (Fig. 1). The PVI 
references canopy reflectance to a soil background line determined 
from regression analysis of the bare soil reflectance measurements. 

The soil background line for the Shiprock sandy loam was 

BRFl - 1.03 BRF2 - 0.05 (r^ - 0.89) [5] 

where BRFl and 2 are bare soil bidirectional reflectance factors In 
channels 1 and 2, respectively. Soil line intersection coordinates, 
determined using the method outlined by Richardson and Wiegand (1977), 
were 


and 


Rgl * -0.024 + 0.515 Rpl + 0.500 Rp2 [6] 


Rg2 = 0.025 + 0.500 Rpl + 0.486 Rp2 [7] | 

from which PVI (Eq. 4) was calculated. 

PVI can also be calculated using radiance. However, since radi- 
ance is proportional to irradiance, a PVI based on radiance will 
change with sun angle (Richardson, 1981)* This can be overcome by 

normalizing the data to a constant sun angle (Jackson et al., 1980). 

Percent cover of the plots was determined using 35 mm color 

slides of the plots (photographed from a vertical position approxi- 
mately 1 m above the plots) projected on a random dot grid. Figure 2 

shows changes In percent cover of the plots during the study. Surface 

soil water content (0 to 4-cm layer) was determined gravimetrical ly on 
samples collected at the time of the spectral measurements. Profile 
soil water content was determined using the neutron attenuation 
method. Meteorological data collected at the site Included dally 
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I, maximum and minimum air temperature, relative humidity, precipitation, | 
2 1 daily solar radiation, and windspeed at 2 m. 

Actual ET was calculated using the water balance equation 
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ET + P - 0 - aSWC [ 8 ] 

where I and p are irrigation and precipitation amounts, respectively, 
0 is drainage, and aSWC is the change in soil water content from the 
previous measurement. Daily ETp was calculated using the Penman 
(1963) equation 

ETp(ly) - (Rn-G) (15.36) (I.O+O.OIM ) (65 - e^) [9] 

Where a is the slope of the saturation vapor pressure-air temperature 
curve, Y is the psychrometric constant, Rn(ly) is net radiation, 
6 (ly) is soil heat flux, es(mb) is the mean saturation vapor 
pressure, ed(mb) is the mean dewpoint temperature, and W(miles) is 
the wind run. 
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The relationship between (defined as ET/ETp) and percent 
cover for the four alfalfa plots is shown in Fig. 3. As expected, 
Kc increased with canopy cover. The data in Fig. 3 include both wet 
and dry soil surface conditions (see Fig. 6 B), as well as advection 
(Kc greater than 1). The linear relationship (r^*0.815) in Fig. 3 
can be thought of as a mean crop coefficient curve for the conditions 
encountered during the study. 

A logical role for remote sensing in estimating Kc is to use 
spectral reflectance to estimate a growth parameter, such as percent 
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cover, and then relate the spectral estimates of growth to Kc using ‘ 
an established relationship such as that of Fig. 3. Although radiance 
or spectral reflectance in properly selected wavebands usually cor- 
relate with crop growth, radiance ratios and vegetation index models 
provide better day-to-day comparisons of spectral data (Richardson and 
Wiegand, 1977). Since soil background reflectance can significantly 
affect reflectance measurements of a crop canopy, we chose to use the 
perpendicular vegetation index (PVI) of Richardson and Wiegand (1977) 
in our analyses. The PVI reduces the effect of the soil background by 
referencing canopy reflectance to a soil background line. A complete 
discussion of vegetation index models can be found in Richardson and 
Wiegand (1977) and Wiegand et al. (1979). 

A simple model can be used to illustrate the behavior of spectral 
reflectance in radiometer channels 1 (0.63 to 0.69 ym) and 2(0.76 to 
0.90 ym), and PVI as a crop canopy grows. Composite bidirectional 
reflectance (vegetation plus soil background) can be approximated by 
the equations 

Rpl = fvRvl + (l-fv)Rgl [10] 
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Rp2 = fvRp2 + (l-fv)Rg2 [11] 
where Rpl and Rp2 are composite BRF's in channels 1 and 2, fv is per- 
cent cover expressed as a fraction, Rvl and Rv2 are vegetation BRF's 
in channels 1 and 2, and Rgl and Rg2 are the corresponding soil back- 
ground BRF's. 

Composite BRF's were calculated for wet and dry soil backgrounds 
(11.5 and 1.0 percent gravimetric water content, respectively) using 
the following values representative of alfalfa and Shiprock sandy 
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loam: Rvl » 0.08, Rv2 « 0,50, Rgl wet * 0.20, Rgl dry * 0.30, Rg2 wet 

» 

» 0.25, and Rg2 dry » 0.35. It was assumed that the soil background 
was sunlit. 

As shown in Fig. 4, composite red reflectance (Rpl) decreases and 
near infrared reflectance (Rp2) increases as percent cover increases. 
The coordinates (Rpl, Rp2) converge to a single point at full cover. 
The area within the triangle in Fig. 4 represents the range of 
possible composite reflectances, assuming that all soil seen by the 
radiometer is sunlit. For a given percent cover, PVI is identical for 
wet and dry soil backgrounds (Table 2). 

In reality, a radiometer will see both shaded and sunlit soil. 
The result will be composite reflectances which are lower than those 
shown in Fig, 4 for the wet and dry soil backgrounds, but values of 
PVI will remain about the same, as shown in Table 3. 

Actual values of PVI were linearly related (r^ = 0.911) to per- 
cent cover of the alfalfa plots (Fig. 5), Th^s compares with an r^ of 
0,861 when a channel 2/1 radiance ratio was compared with percent 
cover. Inclusion of channel 3 (1.55 to 1.75 ^m), a liquid water 
absorption band, in the analyses did not significantly improve the 
results. 

Although Pyi can be used to estimate percent cover and thus a 
mean K<;, it cannot be used as an index to adjust for increased 
evaporation following irrigation since the magnitude of PVI is related 
only to canopy growth. However, the intersection coordinates (Rgl, 
Rg2) on the soil background line, which are used to calculate PVI, are 
related to soil reflectance. Thus, the position of the intersection 
on the soil line may be an indicator of irrigation or rainfall since 
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I wet soil reflectance is less than for a dry soil. As a means of 

j describing the the soil line Intersection, we defined a parameter D as 

( 

D ■ [(Rgl - “o)’' * 


di 

“j where a© is the y-intercept of the soil background line. As shown 

5 ‘ 

in Fig. 1, 0 is the distance along the soil line from the y-intercept 

r ' 

. of the soil line to the intersection coordinates (Rgl, Rg2). 

" i 

*, D and surface soil water content during the study are shown in 

Fig. 6. Irrigations occurred on 15 and 22 July (Table 1). On the 

o! 

I first measurement days following the irrigations, 0 had decreased from 

,0 

the value on the previous day on which spectral measurements were 

,1 

made. The reductions in D ranged from 4.5 to 13.9 percent on 17 July, 
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and from 5,1 to 9.4 percent on 23 July. D then increased as the soil 
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surface dried (Fig. 6), D was sensitive to increases in surface soil 
water content, even at high percent cover. Percent cover on 23 July 
was as high as 90 percent. Differences in D between consecutive 
measurement dates were significantly correlated (r * 0.75^’^) with 
corresponding differences in surface soil water content. Thus, a 
decrease in D from its value on the previous observation date may be 
used as an indicator of irrigation or rainfall so that Kc can be 
adjusted to 1 account for increased evaporation. 

This study has demonstrated a potential use of spectral reflec- 
tance for evaluating the crop coefficient. Once experimental rela- 
tionships between the crop coefficient and stage of development are 
established, repetitive reflectance measurements from aircraft or 
satellites could be used to estimate the stage of development of 
individual fields. Use of reflectance for evaluating crop growth has 
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1 , been denionst rated for a variety of crops (Aase and Siddoway, 1980; ] 

' • I 

2 1 Heilman et al., 1977; Tucker et aU, 1979; Tucker» 1979; Wiegand et 
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3 ! al., 1979). If frequent observations are available, the soil line 

. 1 ' intersection coordinates of PVI could be used as an indicator of when 

s 

5' to adjust the crop coefficient to account for effects of irrigation or 
^ rainfall. 
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Fig. 1, A diagram illustrating the perpendicular vegetation index 

(PVI) of Richardson and Wiegand (1977). A perpendicular 
from the bidirectional reflectance (BRF) coordinates of the 
plant canopy (Rpl, Rp2) Intersects the soil background line 
at coordinates (Rgl> Rg2). The position of Rgl and Rg2 
depends on the reflectance of the soil beneath the crop. D 
is the distance along the soil line from the y«1ntercept of 
the soil line to the PVI intersection coordinates (Rgl, 
Rg2). 

Fig. 2. Changes in percent cover of the alfalfa plots during the 
study. A linearly decreasing amoung of irrigation water was 
applied to the plots with plot 1 receiving the most water* 
and plot 4 the least. 

Fig. 3. Relationship between the crop coefficient (ET/ETp) and 
percent cover of the alfalfa plots. 

Fig. 4. Calculated bidirectional reflectance factors (Rpl and Rp2) 
in radiometer channels I (0,63 to 0.69 pm) and 2(0.76 to 
0.90 pm), respectively, for a wet and a dry soil back- 
ground. Numbers adjacent to the data points denote percent 
cover. 

Fig. 5. Relationship between PVI and percent cover of the alfalfa 
plots. 

Fig. 6. Comparison of D, the distance from the y-intercept of the 

soil background line to intersection coordinates Rgl, Rg2, 
(A) and soil water content in the 0 to 4-cm layer (S). 
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Table 1. Timing and amounts of water applied to the differen- 
tially Irrigated plots of alfalfa, 

I 

I 

?'! 


Date 

Plot No. 

Amount of applied water 
(mm) 

15 July 

1 

98 


2 

70 


3 

50 


4 

38 

22 July 

1 

72 


? 

55 


3 

43 


4 

30 






Table 2 


Simulated values of composite reflectance and PVI for wet 
and dry soil backgrounds (11.5 and UO percent gravimetric 
water content, respectively). 


Percent Cover 


Rpl 

Rp2 

PVI 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

20 

0.18 

0.26 

0.30 

0.38 

0.057 

0.059 

40 

0.15 

0.21 

0.35 

0.41 

0.110 

O.lll 

60 

0.13 

0.17 

0.40 

0.44 

0.163 

0.163 

80 

0.10 

0.12 

0.45 

0.47 

0.216 

0.216 

100 

0.08 

0.08 

0.60 

0.50 

0.267 

0.267 



Table 3. Measured composite reflectances for observed values of 

percent cover comparable to those In Table 2, and associ- 
ated values of PVl. 


Plot No, 

Observed 
percent cover 

Surface soil 
water content 
(% by weight) 

Rpl 

Rp2 

PVl 

3 

19 

1.1 

0.17 

0.26 

0.033 

4 

39 

3.7 

0.13 

0.34 

0.119 

2 

58 

3.9 

0.10 

0.36 

0.154 

1 

80 

3.5 

0,11 

0.43 

0.218 
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Fig. 1. A diagram illustrating the perpendicular vegetation index I 

(PVI) of Richardson and Wiegand (1977). A perpendicular j| 

from the bidirectional reflectance (BRF) coordinates of the I 

plant canopy (Rpl, Rp2) intersects the soil background line 
at coordinates (Rgl, Rg2). The position of Rgl and Rg2 
depends on the reflectance of the soil beneath the crop. D- 
is the distance along the soil line from the y-intercept of 
the soil line to the PVI intersection coordinates (Rgl, Rg2). |1 
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Fig. 3, Relationship between the crop coefficient (ET/ETp) and 
percent cover of the alfalfa plots. 
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Calculated bidirectional reflectance factors (Rpl and Rp2) 
in radiomL:er channels 1 (0.63 to 0.69 and 2(0.76 to 0.90 
pm), respectively, for a wet and a dry soil background. 
Numbers adjacent to the data points denote percent cover. 
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Fig. 6. 


Compapison of D, the distance from the y-intercept of the 
soil background line to intersection coordinate Rqi Ra2 
(A) and soil water content In the 0 to 4-cm layer (B). * * 




